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Abstract. Multi- frequency radio polarimetric observations of the diffuse Galactic synchrotron background enable us 
to study the structure of the diffuse ionized gas via rotation measure maps. However, depolarization will introduce 
artifacts in the resulting rotation measure, most notably in the form of narrow, elongated "depolarization canals". 
We use numerical models of a non-emitting Faraday rotating medium to study the RM distribution needed to 
create depolarization canals by depolarization due to a finite beam width, and to estimate the influence of this 
depolarization mechanism on the determination of RM. We argue that the depolarization canals indeed can be 
caused by beam depolarization, which in turn is a natural consequence when observing a turbulent medium with 
limited resolution. Furthermore, we estimate that beam depolarization can induce an additional error of about 
20% in RM determinations, and considerably less in regions that are not affected by depolarization canals. 
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1. Introduction 

The Galactic magnetic field and turbulence play a ma- 
jor role in structuring the Galactic disk. The recognition 
of the dynamical importance of turbulence in the inter- 
stellar medium (ISM) (e.g. Armstrong et al. H995[) let our 
picture of the ISM evolve from a simple three-phase gas 
to a highly complex medium. Turbulence is believed to 
play a crucial role for molecular cloud and star formation 
(see e.g. Larson lT981l Mac Low & Klessen l2003|) . however, 
its origin is widely debated. On a larger scale turbulence 
leads to chemical mixing (de Avillez & Mac Low 2002 ) and 
contributes to the vertical pressure balance in the Galactic 
disk and to the heating of the ISM (Minter & Balser lT^5|l . 

Magnetic fields are tightly linked to the turbulent na- 
ture of the ISM. Their large-scale components help con- 
fining cosmic rays to the Galaxy, and thus affect chemical 
processes in molecular clouds. They probably contribute 
to generating turbulence in the disk (Sellwood & Balbus 
1999) and definitely affect the vertical disk structure. The 
ratio between mass and magnetic flux in the ISM is an 
essential piece of information necessary to understand the 



role of magnetic fields for molecular cloud and eventually 
star formation (see e.g. Crutcher^)99). 

Density and magnetic field structure in the ISM can be 
probed using multi-frequency radio polarimetry to study 
Faraday rotation 1 in the magneto-ionic medium. Polarized 
extragalactic point sources and pulsars have been used to 
determine rotation measures (RM), but these give only 
measurements at certain lines of sight. Observing the dif- 
fuse Galactic synchrotron background instead yields struc- 
ture in polarized emission on many scales, along many 
contiguous lines of sight. However, depolarization com- 
plicates the interpretation of the resulting RM maps. 
Depolarization occurs if the telescope beam is larger than 
the scale of structure in polarization (beam depolariza- 
tion), along the line of sight if the radiation is emitted and 
Faraday rotated largely in the same medium (depth de- 
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1 Faraday rotation is the rotation of the plane of linear po- 
larization of radiation when it propagates through a magneto- 
ionic medium, due to the birefringence of the medium to circu- 
larly polarized radiation. The rotation of the polarization plane 
through the medium is = RMX 2 , with rotation measure 
RM = 0.81 J B^n e ds, where By is the component of the mag- 
netic field parallel to the line of sight in /iG, n e is the thermal 
electron density in cm -3 and ds is the path length in pc. 
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polarization or internal Faraday dispersion), or if the fre- 
quency bandwidth is so large that the polarization angle 
changes significantly within the band (bandwidth depo- 
larization), see e.g. Burn (1966) or Sokoloff et al. J1998). 
In polarization observations of the Galactic synchrotron 
background, often it is difficult if not impossible to esti- 
mate the separate contributions of beam and depth depo- 
larization. 

Depolarized regions often exhibit canal-like features 
(e.g Duncan et al. 113971 IH^ Gray et al. lTM?! Uyamker et 
al.HHH Haverkorn et al. ESUUUl Gaensler et al. l2"!J0"Tjl . Both 
beam and depth depolarization can cause the canals, but 
require a different underlying RM distribution. The goal 
of this paper is to determine if the observed depolarization 
canals can be due to beam depolarization, and how beam 
depolarization affects the computed RM. To this effect, 
we use numerical models of the warm magneto-ionic ISM 
to compare to observations of the diffuse Galactic syn- 
chrotron background. The models are irradiated with uni- 
form synchrotron background emission, which is Faraday 
rotated in the modeled medium. Then, the polarized ra- 
diation is smoothed to simulate a finite telescope beam, 
where after RM's are computed from the smoothed emis- 
sion maps. The structure in polarized intensity and RM 
emerging from the smoothed modeled radiation field is 
then compared to the observations. This allows us to sim- 
ulate purely the effect that beam depolarization has on 
polarimetric observations of the ionized ISM. 

This work was motivated by the observations briefly 
presented in Sect. [3 In Sect. we describe the numeri- 
cal models, the addition of the synchrotron background 
and the simulation of beam depolarization. Section0]dis- 
cusses the signature of beam depolarization in the maps 
of polarized intensity, while in Sect. El the effect of beam 
depolarization on RM is studied. 

2. The observations 

The observations were made with the Westerbork 
Synthesis Radio Telescope (WSRT) in five frequency 
bands centered on 341, 349, 355, 365, and 375 MHz, each 
with a bandwidth of 5 MHz. The field discussed here is 
centered at (l 7 b) — (161°, 16°) and about 50 square de- 
grees in size, and has a resolution of about 4'. All four 
Stokes parameters /, Q, U and V are measured, from 
which polarized intensity P and polarization angle <f> were 
determined as 

p = VQ 2 + u 2 , (i) 

<f> = — arctan — . (2) 

The instrumental linear polarization is Si 1%. For a de- 
tailed description of this field see Haverkorn et al. ( 2003a). 

No signal was detected in total intensity /. As the in- 
terferometer is only sensitive to scales smaller than ap- 
proximately a degree, the non-detection of I means that 
/ does not exhibit any structure on scales smaller than 
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Fig. 1. Observed polarized intensity at 349 MHz in a field 
centered at (l,b) = (161°, 16°). High intensity is shown 
white, up to a maximum of ~ 13 K. 

~ 1°, so that the observed total synchrotron power must 
be very uniform. On the other hand, abundant structure 
in linearly polarized emission was detected, both in P (see 
Fig.^) and in <f>. This indicates that the small-scale struc- 
ture in polarization is caused by Faraday rotation and 
depolarization. Ubiquitously present in the field are one- 
beam wide "canals" of beam depolarization, which we will 
discuss in Sect.0] Missing large-scale structure in P prob- 
ably does not play a large role in this field (Haverkorn et 
al. 2003a, 2004b). 

RMs were computed from linear fits of 4> against A 2 . 
For about 70% of the beams with an average polarized in- 
tensity (P) > 5 S/N, reliable RM's (with reduced x 2 < 2) 
could be determined. This is about 28% of all data. Fig. [21 
shows a grey scale map of RM, where in each beam a 
RM value is plotted, regardless of the quality of fit of the 
0(A 2 )-relation. In general, the positions showing anoma- 
lous RM with rapidly changing magnitude and sign are 
the positions where noise dominates, so that reliable RM 
determination is not possible, see Sect. 01 

Contributions to the RM are made along the line of 
sight at positions with a non-zero electron density and 
magnetic field. This is the case in the warm and hot ISM 
components. The filling factor in the warm ISM is approx- 
imately 20% (Reynolds HMD, while the elect ron density 
in the hot ISM is so low that its RM contribution can be 
neglected for our purposes. Furthermore, polarized emis- 
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Fig. 2. Map of RM in the observed field. The grey scale 
runs from —15 to 15 rad m~ 2 . The map shows an RM 
value at every position, regardless of the quality of the 
linear fit of <j> against A 2 . 

sion is increasingly depolarized when emitted at greater 
distances. Haverkorn et al. <|2004a[) estimate a distance of 
i~ 600 pc between observer and the "polarization hori- 
zon" , defined as the point beyond which less than about 
a third of the polarized emission is not depolarized while 
propagating through the medium. 

3. Numerical models 

In order to estimate the effect of beam depolarization on 
the RM determinations, we generate artificial maps of po- 
larized intensity at five wavelengths from numerical mod- 
els of a turbulent magnetized ISM. 

3.1. Domain considerations 

Ideally, the numerical models should resemble the obser- 
vational domain (i.e., the ratio between the length of the 
line of sight and the extent in the plane of sky) while 
retaining a resolution in the plane of sky comparable to 
the observations. At an estimated polarization horizon of 
600 pc, the spatial extent of the 7° x 9° field of obser- 
vation is about 75 pc. The observational resolution of 4' 
would correspond to approximately 120 x 120 resolution 
elements in the plane of sky. 



However, for the numerical models, this estimated re- 
quirement is too optimistic for two reasons. First, in order 
to investigate beam depolarization, we need substructure 
within one beam (i.e. resolution element). Second, any 
magnetohydrodynamics (MHD) code suffers from numer- 
ical dissipation at the smallest scales, thus introducing a 
lower limit for physically reasonable structures in the sim- 
ulations (between 2 to 4 resolution elements in our case). 
Both requirements taken together lead to a desirable res- 
olution of ~ 500 x 500 elements in the plane of sky for the 
simulations, resulting in a domain size of 500 x 500 x 4000 
elements. This is beyond current computer capacities. 

Thus we decided to use a cubic simulation domain 
(at 512 3 elements) and to construct a bar-like domain by 
stacking this cube for 8 times, each time rotating it and 
shifting it by one fifth of the domain size. Thus we get dif- 
ferent realizations of the turbulent structures, mimicking 
a turbulent medium over the whole length of the line of 
sight. 

3.2. Model description 

Our "observational" domain is derived from a periodic- 
box-simulation of the turbulent magnetized ISM by Li 
et al. <j2T)07)l HSS|| and Heitsch et al. ffiWty . employ- 
ing ZEUS-MP, the massively parallel version of ZEUS-3D 
(Norman I2000|l . The electron density and the magnetic 
field are rescaled to the values for the warm ISM, from 
which five models with different path lengths through the 
medium were constructed, as shown in Tabled 

The density and the magnetic field are initially uni- 
form, and will be perturbed by a fixed velocity field. We 
generate the velocity field as described by Mac Low (JT999 ) , 
as follows. Each wave number k is given an amplitude in 
Fourier space, which is randomly drawn from a Gaussian 
distribution around unity, and a random phase. Fourier 
transforming to velocity space yields velocity components 
in three dimensions for each position. This is the velocity 
field with which the model is perturbed at each timestep, 
to simulate turbulent driving. The velocity components 
are multiplied by a velocity amplitude at each timestep, 
chosen such that the input energy stays constant with 
time. Perturbations are only induced at the lowest spa- 
tial wavenumbers 1 < k < 2. In the plane of sky, this 
corresponds to the largest scales possible in the domain. 
The driving mechanism is meant to mimic energy input on 
the largest scales by an unspecified physical process such 
as supernova shock fronts or Galactic shear. 

We use a snapshot of the simulation at a time when 
the full turbulent cascade has developed and a steady state 
between energy input on the largest scales and numerical 
dissipation on the smallest scales has been reached. This 
steady state corresponds to a Mach number of M. w 10. 
This is more than the values usually assumed for the warm 
ISM (namely M. « 1 at a sound speed of ks 6 km s^ 1 ), 
thus rendering the turbulence more compressive than in- 
tended for this work. 
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However, the turbulent cascade represented by a 
Kolmogorov slope of —5/3 does not change substantially 
(see e.g. Cho et al.[2003), so that we can rescale the densi- 
ties and held strengths to values appropriate for the warm 
ISM. Fluctuations in electron density are allowed up to 
an upper limit of 10% of the mean density (assumed to 
be 0.1 cm -3 ), which corresponds to the upper limit to 
structure in emission measure (EM) in the held of obser- 
vation given by high-resolution Ha observations (Madsen, 
private communication). The held variations are on the 
order of the mean held of 2.5 /j,G. 

We assume a filling factor of ionized gas of 100%, and a 
sharp cutoff of Faraday rotation of the polarized emission 
at 600 pc. As the RM is an integral over the line of sight, 
a smaller filling factor combined with a longer path length 
would give the same results. 

The MHD model does not contain synchrotron emis- 
sion within the medium itself. This is obviously not a real- 
istic situation, but allows us to investigate the influence of 
beam depolarization without depth depolarization, as the 
latter only occurs in a medium that both Faraday rotates 
and emits. Including synchrotron emission in the medium 
would result in additional depolarization and larger devi- 
ations from the linear </>(A 2 ) relation. 
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Table 1. Parameters of the numerical models used: N x , 
N y and N z are the model dimensions in resolution ele- 
ments and Li os is the length of the line of sight in parsecs. 
The initial magnetic field B\\ is oriented along the line of 
sight z and n e is the initial uniform electron density. ARM 
is the width a of the Gaussian distribution of RM, scaled 
via L. 



domain. Stokes Qo and Uq are derived from those maps 
via 

4> = RM A 2 , (4) 

U = sin(2<£), (5) 

Qo = cos(2(t>). (6) 

Note that we assume 100% uniform background polariza- 
tion, P = vQo + ^o = 1> an d a constant background 
polarization angle of </>o = 0. The wavelength A is given in 
meters. 

To mimic the limited telescope resolution, we smooth 
the Stokes Qo an d Uq maps with a Gaussian beam of width 
a. This results in maps of Q and U which are a times 
oversampled and which now can be used like observational 
maps. A noise level of 0.05 Q, U rms is added to the Stokes 
parameters, which corresponds to the noise in the high 
signal-to-noise regions in the observations. P and <fr are 
calculated according to Eqs. {""J) and J""}. From the angle 
maps corresponding to the five observational wavelengths 
we finally can determine RMji t as 

<^ = RM /lt A 2 . (7) 

We apply the same criteria to select reliable fits as in the 
observations (see end of jJ""J. 

Fig. |31 gives an example for model BA at smoothing 
a = 4. P shows elongated regions of depolarization, not 
directly related to features in the RM map. However, the 
regions of low polarized intensity P occur at locations of 
large gradients in RM. These structures reappear in the 
fitted RM, since regions of low P are noise-dominated 
and thus do not allow an accurate ht. The other mod- 
els show similar characteristics: elongated depolarized re- 
gions which are one beam wide are present in all models, 
and the RM derived from smoothed data is patchy and 
shows sub-beam scale anomalies. However, a larger width 
of the RM distribution ARM induces more depolarization, 
so that larger depolarized regions exist in those models. 
For further distinction between models, see Sect. 03 



3.3. Simulating beam depolarization 

The goal is to produce maps of Stokes Q and U from the 
models, which then can be treated like observational data. 
From the bars we construct the numerical RM in the plane 
of sky as 

L Nz 

RM^ = 0.81-^ 5>eS||k^= <■.! ■ !'--- V,}. (3) 

z k=i 

in units of rad m~ 2 . The line-of-sight length Li os is given in 
parsecs, the electron density n e in cm -3 , and the magnetic 
field component along the line-of-sight Bn in /iG. N X ,N Z 
are the number of cells in x and z direction, respectively, 
and N y — N x . For simplicity, we integrate along parallel 
lines of sight instead of an integration over a cone-like 



4. Depolarization canals 

All high-resolution observations of the polarized intensity 
of the diffuse synchrotron background exhibit elongated, 
one-beam-wide "canals" of very low or zero P, where the 
radiation is depolarized (e.g Duncan et al. 119971 [1999 
Gray et al. 11*3531 Uyam ker et a l. lT3331 Gaensler et al. l*"~"7~Tl 
Haverkorn et al. I2003al I2003b|) . The polarization angle al- 
ways changes by 90° across such a canal (Haverkorn et al. 
2000 ) . Two possible explanations exist for the depolariza- 
tion canals (see e.g. Burn l""5"*""l Sokoloff et al. [13398): 

1. beam depolarization: if the gradient in RM VRM = 
(2n + 1)tt/2 within one beam, the radiation within 
the beam is depolarized. This mechanism would ex- 
plain the fact that canals do not shift in position with 
frequency. However, it requires a medium with steep 
gradients in RM and/or intrinsic angle 4>q within one 
beam, elongated perpendicularly to the gradient. 
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Fig. 3. Polarized intensity P (4xl(T 3 < P < 0.97), gradient of RM, |VRM| (lxltr 3 < \VRM\ < 4.9 rad m -2 pixel -1 ) 
and numerical and fitted RM (scaled identically between —6.8 < RM < 22.32 rad m -2 ) for model B4 at smoothing 
cr = 4, four times oversampled. Note that depolarized (black) canals appear at the locations of large (white) gradients. 
These regions reoccur in anomalous RM fa , leading to fringes in the fitted RM map. 



depth depolarization: a uniform synchrotron emitting 
magneto-ionic medium can cause complete depolariza- 
tion for |RM| = n7r. However, in this explanation the 
depolarization canals shift in position with frequency, 
and a uniform medium (across a major part of the path 
length) is needed for complete depolarization. 



Haverkorn et al. |2004b) argue that in their observations 
(one of which is used here), beam depolarization due to 
RM gradients is the dominant mechanism causing depolar- 
ization canals. This does not in any way mean that depth 
depolarization is not important in the rest of the medium. 
For observations made at 350 MHz, only a modest gra- 
dient in RM across a beam (~ 2.1 rad m~ 2 ) can create 
a canal, whereas much larger gradients (~ 34 rad m~ 2 ) 
are needed to form a canal in 1.4 GHz observations. For 
an interpretation of canals at 1.4 GHz in terms of depth 
depolarization, see Shukurov & Berkhuijsen ( 2003). 



In the numerical models discussed in Sect. depo- 
larization canals are visible as well (see top left map in 
Fig. [3J , which are caused by beam depolarization. 

At first sight, there is no obvious connection between 
the locations of the canals and RM (left hand maps in 
Fig. 0J- However, the magnitude of the gradient in RM 
(top right map in Fig. |3J) is obviously correlated to the po- 
sition of the canals. This is demonstrated more quantita- 
tively in Fig. 0] which shows a scatter plot of P against the 
magnitude of the gradient in RM with smoothing a = 4 
for each independent "beam". The solid line denotes the 
theoretically expected depolarization in a beam due to a 
resolved gradient dRM/dr in RM (Sokoloff et al. H355) l 



P 

— =exp 

■M) 



_L_ /gMV 4 
'in 2 \ dr ) 
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where Pg is the polarized intensity for constant RM within 
the beam. The polarized intensity decreases almost to zero 
where the gradient in RM is such that A(f> ~ w/2. 

We selected "canal-like beams" by comparing P in a 
certain point to the average P of two diametrically op- 
posed neighboring beams. If P in the central beam was 
less than 20% of the averaged P of the neighboring beams, 
the beam was considered a canal. Note that this selection 
criterion does not use any information about the length 
of the canals. Fig. shows the distribution of the gradi- 
ent in RM only at canals, defined according to the above 
process. The five plots show canals defined in five wave- 
lengths, and the solid and dotted lines represent smooth- 
ings cr = 4 and 8, respectively. The dashed lines denote 
the value of VRM at the location where Acf> = —ir/2, w/2. 
Clearly, the gradient in RM is peaked at the value where 
A<p = tt/2. The decline in number of canal pixels with fre- 
quency can be explained by the fact that the width of the 
polarization angle distribution decreases with increasing 
frequency for a given RM distribution. 

So the gradient of RM in the canals has the right mag- 
nitude to cause beam depolarization (Fig. |SJ. Since the 
canals exist over many beams, the most likely direction 
of this gradient is perpendicular to the direction of the 
canals. In Fig. we show P in grey scale, with vectors 
of the RM gradient superimposed, demonstrating that in- 
deed the gradients in RM are directed perpendicular to 
the canals. 

Detailed analysis of shape, direction and origin of RM 
gradients and the implications on the medium will be dis- 
cussed in a forthcoming paper. However, we give a simple 
estimate here of the required changes in magnetic field 
strength B and/or electron density n e to produce large 
enough gradients. For A = 0.80 m, ARM = 2.45 rad m~ 2 
in order to get an angle change of A<j) — n/2 within a 
beam. If the structures in VRM are filaments at an aver- 
age distance of 300 pc, their width would be 0.4 pc. The 
necessary ARM requires an enhancement of n e and B of a 
factor 3 to 4, for ISM parameters as in e.g. model B2 (see 
Table ^ . On the other hand, if the structures are sheet- 
like, and have an extent along the line of sight of a few 
parsecs, then an increase in electron density and magnetic 
field of 50% is sufficient, or an increase in either n e or 
B of 100%. However, the unpublished Ha measurements 
discussed above (Madsen, private communication) give an 
upper limit to small-scale structure in emission measure 
EM = J rig dl « 0.3 cm~ 6 pc, which corresponds to a max- 
imum An e w 0.02 cm -3 . Therefore, the gradients in RM 
are probably mainly caused by magnetic field changes, and 
could originate from shocks, local narrow magnetic field 
reversals or other sheet-like structures in density and/or 
magnetic field. 

Note that the canals in the numerical models are much 
shorter and less ordered than the ones in the observations. 
This is because in the numerical models we can introduce 
dynamical structure self-consistently only up to the do- 
main size. Ideally, we would like to use subframes from 
the numerical models to be able to include larger scale 
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Fig. 4. Absolute value of the gradient in original RM in 
the simulations |Vi?Af| against polarized intensity P of 
the smoothed maps with smoothing a — 4 for each pixel. 
The curve corresponds to Eq. [SJ 

structure in the models, so that modeling of canals across 
a significant part of the field would be possible. However, 
the available resolution is too low for this to be feasible. 

Elongated narrow gradients in magnetic field and/or 
density have also been seen in other numerical simulations 
(e.g. Schekochihin et al. l2Tjo"2l Cho et al. l2"U0ljl . 

Concluding, in the numerical models discussed here 
narrow and elongated gradients in RM exist of large 
enough magnitude to cause beam depolarization in the 
modeled medium. Taking into account other observational 
evidence (Haverkorn et al. I2004bj) . this leads us to con- 
clude that the narrow and deep depolarization canals as 
seen in the observations are most likely caused by beam 
depolarization due to RM gradients. 

5. Reliability of RM determinations 

The beam depolarization raises the issue of how reliable 
the RM determinations are. Fig. quantifies the effect 
of limited resolution on the values of the computed RM. 
We determine the normalized mean quadratic difference 
between fitted and original smoothed RM as 



*RM=(» 



RM /it ) 2 > 



(RM ) 



1/2 



(9) 



where RM is the original RM and RM CT is RM smoothed 
over Gaussian beams of width a. We have chosen for a 
<5RM defined by Eq. © for the following reason. The goal 
of the equation is to compare the fitted RM to the "real" 
RM in the modeled medium. On scales below the beam 
size, the fitted RM cannot show any structure that repre- 
sents structure in the medium by definition. Therefore, if 
RMfu would be compared to the RM as directly derived 
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Fig. 5. Distribution of difference in RM across a canal, ARM. From left to right, canals are selected in wavelengths 
of 88cm, 86cm, 84cm, 83cm and 80cm, respectively. The solid line denotes canals in P maps smoothed by a — 4, the 
dotted lines are for a smoothing a = 8. The dashed lines give the positions where ARMA 2 = — it/2, n/2. 




Fig. 6. Part of the polarized intensity map at smoothing 
a = 4 in grey scale, with strength and direction of the 
gradient in RM superimposed as vectors. 
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Fig. 7. Deviation of RM, <5RM, calculated from smoothed 
Q and U maps from original RM, relative to original RM, 
against smoothing a (see Eq. EJ . Different symbols denote 
results for models Bl to B5 with a given width of the 
original RM distribution. Only independent beam posi- 
tions with (P > 5S/N) A (x 2 < 2) are used to determine 
SUM. 



in the model, RMfu — RM would give large values which 
are not the result of the smoothing, but only due to the 
fact that the beam size is incorporated in the determina- 
tion of RM fn, but not in that of RM. Therefore, we chose 
to compare RMj, f to an RM which is smoothed over a 
Gaussian beam, assuming that the smoothed RM approx- 
imates the physical RM in the model if viewed only on 
scales larger than the beam size. In addition, the subse- 
quent normalizing to the distribution width of the original 
unsmoothed RM rules out any spurious effects on the dif- 
ference estimates due to averaging. In this way, SRM is 
a measure of how accurately the RM computed from Q 
and U data with a finite beam width approximates the 
physical RM in the medium. 

Fig. [7] shows the deviations <5RM for all independent 
beams satisfying the criteria for good fits (see end of 
For a width of the RM distribution ARM « 20 rad m -2 



and a smoothing of a = 4 (values which are typical for the 
observations by Haverkorn et al. 2003a, 2003b) we would 
expect from Fig. an error of ~ 20% in the fitted RM. 
However, as we will see below, this error estimate is an 
upper limit. 

For ARM small enough not to lead to sign changes of 
Q and U within a beam, the differences SRM saturate at 
a 5%-level with increasing beam width, since variations of 
RM within a beam are limited (Fig.0for ARM = 3). This 
behavior changes profoundly if Stokes Q and U reverse 
sign within one beam. This is equivalent to a change in 
polarization angle of 90° which, as discussed above, leads 
to depolarization canals. The effect on SRM is twofold: 

1. Low P means small S/N, thus not well defined angles, 
and thus large errors in the fits. Therefore, those fits 
are excluded by the selection criterion. 
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2. Averaging introduces slight non-linearities in the an- 
gles around depolarized regions, i.e. regions where 
Stokes Q and U reverse sign within one beam. These 
non-linearities lead to an overshoot of the RM (see 
Fig. [SJ close to the depolarized region. 

With increasing ARM, Stokes Q and U are more 
likely to change sign within a beam. Beam averaging 
then leads to much severer depolarization (see Fig. [SJ. 
Depolarized regions are accompanied by systematic over- 
shoots in RM f u (Fig. [SJ top panel) , which stem from the 
non-linearities in Eq. @ arising when Stokes Q changes 
sign. (Note that in the bottom panel of Fig. [SJ the polar- 
ization angles derived from beam-averaged Stokes Q and 
U are not ordered with respect to A any more in a canal, 
thus leading to ill-defined fits.) The overshoots of extreme 
positive and negative RM are also visible in the observa- 
tional data (see Fig. El, indicating that the depolarization 
canals in the data are also caused by beam depolarization 
as opposed to depth depolarization. 

Not all RM gradients introduce these spurious results. 
Far away from the gradient, the fits retrieve the original 
RM. However, close to the gradient, multiple oscillations 
of Qq and Uq within the gradient lead to Q and U close 
to zero over an extended region. 

For one isolated RM-gradient as shown in Fig.[SJi (see 
figure caption), the differences SKM selected for reliable 
fits behave with the change in RM across the canal ARM 
as shown in Fig. The first spike corresponds to a Oth 
order canal with an angle difference of ir/2. Higher-order 
canals have angle differences of (2n + 1) 7r/2, as discussed 
above. The RM field from the numerical models is a com- 
bination of all Ai?Ms given in Fig. [5J 

Note that the fitted RM/it deviates from the origi- 
nal RM most for those RM gradients where canals occur. 
In between the spikes, at positions where canals are not 
near, the difference <5RM drops nearly to zero, i.e. a fi- 
nite beam width has almost no effect on the RM deter- 
mination. Therefore, 5RM is dominated by values around 
canals where P > 5 S/N (needed to be selected as reliable 
RM). 

From Fig. El we conclude that the difference estimates 
5RM for ARM > 3 rad m -2 as given in Fig. are indeed 
too pessimistic, since they are dominated by the regions 
around canals, that are influenced by the canals but still 
fall within the definition of "reliable RMs" . For smooth re- 
gions, or regions not affected by canals, we expect the RMs 
determined from the observations to be accurate within 
5% of ARM, corresponding to an error of 0.15 rad m~ 2 
for the observations mentioned here. The overall error will 
increase with increasing number of canals, which is why 
Fig. generally shows larger SRM than Fig. El 

6. Conclusions 

Multi-frequency radio observations of the Galactic syn- 
chrotron background allow us to probe the structure of the 
magneto-ionized component of the ISM from RM maps. 
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Fig. 8. Details of an idealized canal caused by sign changes 
in Stokes Q (see text), (a): RM against spatial coordinate 
x. The original RM (thin line) has a gradient between 
the two vertical dotted lines, leading to a smoothed RM S 
(thick solid line) and a fitted RM/jj (thick dashed line) 
showing an overshoot, as in the observations. Smoothing 
("beam") width corresponds to Ax = 4. ARM as used in 
the text is the change of RM over the gradient, in this case 
ARM = 6.5. (b)-(d): Stokes Q, U and polarized intensity 
P for all wavelengths before smoothing (thin lines) and 
after smoothing (thick lines): Strong oscillations of the 
unsmoothed Q and U within the gradient region of RM 
(see Eqs. & [HI) lead to cancellation of smoothed Q and 
U within that region. Correspondingly, P drops to nearly 
0. (e): Polarization angles for all wavelengths (as denoted 
in legend of Q). Thin lines show the angles derived from 
the original, unsmoothed RM, thick lines denote the angles 
derived from the smoothed Q and U. Note that within the 
gradient region, there is no ordering of <p with A. 

However, the quality of the RM determinations depends 
strongly on the degree of local depolarization due to in- 
trinsic structure in the ISM, and thus in the RM. We in- 
vestigated the effect of beam depolarization on the quality 
of the fitted RMs with the help of numerical models. We 
find: 
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Fig. 9. 6RM (see Eq. EJ) for an isolated RM-gradient as 
shown in Fig.lBJa) against ARM, the change of RM across 
the gradient. The gradient width is 4 pixels while the beam 
width varies between 4 and 10 pixels. Fitted RMs are 
derived in the same way as for the models, selected for 
(P > 5S/N) A (x 2 < 2) and independent beams. Note 
the periodicity T RM w 4.4 in <5RM, corresponding to an 
angle change of ir. The first spike occurs at an angle dif- 
ference of 7r/2, the canal of 0th order. 

1. The models show depolarization canals due to beam 
depolarization across sharp gradients in RM. This sug- 
gests that narrow gradients of a few rad m -2 can 
be a common feature in a magneto-ionized medium. 
Combined with other observational evidence, we there- 
fore conclude that the depolarization canals seen in 
our 350 MHz observations are caused by beam depo- 
larization due to RM gradients. While the width of 
the canals is an observational artifact, their presence 
and length mirror underlying structure in the Faraday- 
rotating medium. 

2. RM determinations close to canals can be incorrect 
by several rad m~ 2 , depending on the absolute RM 
difference between both edges of the canal. This is a 
consequence of a non-linearity in 4>(X 2 ) introduced by 
finite beam width. 

3. RM determinations in smooth regions or in regions 
separated from canals by at least a beam are accurate 
to 0.05ARM, where ARM is the width of the RM dis- 
tribution, or SRM = 0.15 rad m -2 for the observations 
mentioned here. 

In future research, improved models will enable us to 
investigate what structures in the ISM actually cause the 
elongated gradients of RM, and will serve as a calibration 
what structure information we can extract from observa- 
tional maps. 

Acknowledgements. The authors wish to thank Ellen Zweibel, 
Peter Katgert, Mordecai-Mark Mac Low and Enrique Vazquez- 
Semadeni for enlightening and invigorating discussions. We are 



grateful to P. S. Li and M. L. Norman for making the 512 
available, and thank Wolfgang Reich for useful and critical 
comments. Computations were performed on the SGI Origin 
2000 of the National Center for Supercomputing Applications 
(NCSA) at Urbana-Champain. The Westerbork Synthesis 
Radio Telescope is operated by the Netherlands Foundation 
for Research in Astronomy (ASTRON) with financial sup- 
port from the Netherlands Organization for Scientific Research 
(NWO). The project originated during the Astrophysical 
Turbulence Program at the KITP at University of Californa, 
Santa Barbara in 2000 (NSF grant PHY94-07194). M.H. is 
supported by NWO grant 614-21-006. F.H. is supported by 
a Feodor-Lynen fellowship of the Alexander-von-Humboldt 
Foundation and by NSF grant AST-0098701. 

References 

Armstrong, J. W., Rickett, B. J., & Spangler, S. R. 1995, ApJ, 
443, 209 

de Avillez, M. A., & Mac Low, M.-M. 2002, ApJ, 581, 1047 
Burn, B. J. 1966, MNRAS, 133, 67 

Cho, J., Lazarian, A., & Vishniac, E. T. 2002, ApJ, 564, 291 
Cho, J., Lazarian, A., & Vishniac, E. T. 2003, ApJ, 595, 812 
Crutcher, R. M. 1999, ApJ, 520, 706 

Duncan, A. R., Haynes, R. F., Jones, K. L., & Stewart, R. T. 

1997, MNRAS, 291, 279 
Duncan, A. R., Reich, P., Reich, W., & Fiirst, E. 1999, A&A, 

350, 447 

Gaensler, B. M., Dickey, J. M., McClure-Griffiths, N. M., 
Green, A. J., Wieringa, M. H., & Haynes, R. F. 2001, ApJ, 
549, 959 

Gray, A. D., Landecker, T. L., Dewdney, P. E., Taylor, A. R., 
Willis, A. G., Normandeau, M. 1999, ApJ, 514, 221 

Haverkorn, M., Katgert, P., & de Bruyn, A. G. 2000, A&A, 
356, LI 

Haverkorn, M., Katgert, P., & de Bruyn, A. G. 2003a, A&A, 

403, 1031 

Haverkorn, M., Katgert, P., & de Bruyn, A. G. 2003b, A&A, 

404, 233 

Haverkorn, M., Katgert, P., & de Bruyn, A. G. 2004a, submit- 
ted to A&A 

Haverkorn, M., Katgert, P., & de Bruyn, A. G. 2004b, submit- 
ted to A&A 

Heitsch, F., Zweibel, E. G., Mac Low, M.-M., Li, P. S., Norman, 

M. L. 2001, ApJ, 561, 800 
Larson, R. B. 1981, MNRAS, 194, 809 

Li, P. S., Norman, M. L., Heitsch, F., & Mac Low, M.-M. 2000, 

BAAS, 197, 502 
Li, P. S., Norman, M. L., Mac Low, M.-M., & Heitsch, F. 2004, 

ApJ, accepted 
Mac Low, M.-M. 1999, ApJ, 524, 169 

Ma c Low, M.-M., fc Klessen, R. 2003, ARAA, in press, 

astro-ph/0301093 
Minter, A., & Balser, D. S. 1998, LNP, 506, 543 
Norman, M. L. 2000, Revista Mexicana de Astronomia y 

Astrofisica Conference Series, 9, 66 
Reynolds, R. J. 1991, ApJ, 372, 17L 

Schekochihin, A. A., Maron, J. L., Cowley, S. C, & 

McWilliams, J. C. 2002, ApJ, 576, 806 
Sellwood, J. A., & Balbus, S. A. 1999, ApJ, 511, 660 
Shukurov, A., & Berkhuijsen, E. M. 2003, MNRAS, 342, 496 
Sokoloff, D. D., Bykov, A. A., Shukurov, A., Berkhuijsen, E. 

M., Beck, R., & Poezd, A. D. 1998, MNRAS, 299, 189 
Uyamker, B., Fiirst, E., Reich, W., Reich, P., & Wielebinski, 

R. 1999, A&AS, 138, 31 



